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Abstract 



Entanglement between two quantum systems is a resource in quan- 
tum information, but dissipation usually destroys it. In this article 
■ we consider two qubits without direct interaction and we show that, 

even in cases where the open system dynamics destroys any initial en- 
tanglement, the mere monitoring of the environment can preserve or 
create the entanglement, by filtering the state of the qubits. While 
the systems we study are very simple, we can show examples with 
entanglement protection or entanglement birth, death, rebirth due to 
monitoring. 

Keywords: Entanglement; Concurrence; Dissipative dynamics; Con- 
tinuous observation; A priori state; A posteriori state. 

1 Introduction 

Entanglement is an intrinsically quantum type of correlation among quan- 
tum systems which is of fundamental importance in quantum information 
[1]. The behaviour of entanglement under dissipative dynamics has been 
studied extensively [3l[2], either to find means to protect entanglement 
against decoherence, either to understand how to use a dissipative dynamics 
to create entanglement. Usually dissipation tends to destroy entanglement, 
at least when the two quantum systems do not interact directly. Sometimes 
this disentaglement can be completed even in a finite time [HEllS] and this 
phenomenon has been called entanglement sudden death (ESD). However, 
dissipation can create entanglement too; this happens when the two parties 
interact with a common bath [BHSIEIE] ) even if they do not interact directly, 
and we can have entanglement birth, death, rebirth. Entanglement can be 
preserved or generated also by controlling the composite system by means 
of measurement based feedback P-fTT]. 

Preservation of entanglement can be obtained also by pure monitoring 
of the system |12H14j . that is by an indirect measurement on the system 
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which acquires information thanks to the observation of its environment, 
but which does not perturb the system. Quantum trajectory theory ahows 
for describing a continuous monitoring |15pi6j and in such a theory we have 
to distinguish between the a posteriori state, the conditional state given 
the observed output, and the a priori state, the mean state, satisfying a 
master equation. It is possible that the a posteriori states are entangled, 
while the a priori state is not. By using the concurrence [T7] as a measure 
of entanglement it has been shown that the pure monitoring can slow down 
the decay of the entanglement [12] . 

The aim of our paper is indeed to study the effect of monitoring on the 
a posteriori entanglement when the a priori dynamics washes out any initial 
entanglement. More precisely, we consider the case of the open dynamics 
of two qubits in the Markovian regime and we model their global evolu- 
tion by a Hudson-Parthasarathy equation. This approach allows to clearly 
characterize the Markovian evolutions representing two qubits which do in- 
teract or do not interact, directly or through a common bath. Section [2] 
is devoted to the HP evolutions and to such a characterization; we recall 
also how to introduce measurements continuous in time and how to get the 
corresponding stochastic Schrodinger equation (SSE) and stochastic master 
equation, which are the starting points to study the dynamical behaviour 
of the monitored system and of its entanglement. In Section [3] we consider 
the case of no direct or indirect interaction between two qubits. When only 
local detection operators are involved, we show that, by pure monitoring, 
the decay of entanglement can be slowed down and, in special cases, even 
stopped independently of the qubit initial state (entanglement protection). 
In cases with non local detection operators, we show that, now depending 
on the qubit initial state, entanglement can even be created by pure moni- 
toring (entanglement generation) . In Section H] we study a case of indirect 
interaction between the two qubits through a common bath. We show that, 
even if the a priori dynamics completely destroys any entanglement, a proper 
monitoring scheme can maximally entangle any initial qubit state. 

1.1 Two qubits 

We consider two qubits; for each qubit we denote by |1) the up state and 
by |0) the down state. By ax, Cy, Uz we denote the Pauli matrices. In 
^ = the canonical basis (or computational basis) [1] is 

|wi) = |ll), |U2) = |10), |U3) = |01), |U4) = |00), (1) 

and the Bell basis [18] is 

|/3o) = -^(|00) + |ll)), lA) =(7,^a|/3o), i = l,2,3. (2) 

The set of statistical operators is S(!K) and the one of linear operators 
is £(^). A local operator is a linear operator which acts non trivially only 
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on one of the factors of (E> C^, i.e. it has the form Ti or Ti A with 
A £ £(C^). The two qubits are independent if their state is a product state 
p = Pi ® P2- The separable states [19] are the statistical operators which 
admit a convex decomposition into product states, so that the correlation 
between the two qubits has a classical explanation; the other statistical 
operators are said to be entangled. The maximally entangled states are the 
pure states which, by partial trace on one of the two factors, reduce to 
maximally chaotic states, that is 1/2. The projection on one of the Bell 
vectors ([2]) is a maximally entangled state. 

1.2 Concurrence 

A very useful measure of entanglement is the concurrence, introduced by 
Wootters ^\ . Let us consider a generic vector (p £'K and expand it on the 
canonical basis ([1]) 

if = (/Jiilll) + vPiollO) + (/^oilOl) + (/^oolOO). (3) 

Let T be the complex conjugation of the coefficients in the canonical basis: 

T(/p = ^|ll)+^|10)+m|01)+^|00). (4) 

Let us define 

■= (T9?ky ® ayif) = 2 {pio^m - ipiipoo) , := \x^\ . (5) 

When \\(p\\ = 1, Cp is the concurrence of the pure state ip. In general, if p 
is not normalized and = jj^, then 

= (6) 

\m\ 

Note that Cf^^ = 1 and Cu^ =0. 

If /? is a generic statistical operator, the concurrence is defined by 

Cp :=inf^p,C^^, (7) 

i 

where the infimum is taken over all decompositions of p in pure states, p = 
see for instance [3] p. 231. We have < Cp < 1, Vp S S(^), 
with Cp = if and only if p is separable and Cp = 1 if and only if p is 
maximally entangled. 

A subclass of states, for which it is easy to compute the concurrence, is 
the one of the "X" states [5l|3]: in the canonical basis, an X state has non 
vanishing matrix elements only in the two main diagonals. The projection 
on a Bell vector is an X state. For any X state p, by setting pij = {ui\puj), 
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we have pjj > 0, pij = pji, J2j=iPjj = 1) PiiPu > \pu\ , P22P33 > IP23I ; 
moreover, the concurrence is given by [5] 

Cp = 2 max {0, Ci, C2}, (8a) 

Cl = \p23\ - ^/PllPM, C2 = \Pia\ - \J P22P33- (8b) 

FinaUy, let A and B be hnear operators on C^. In studying the dynamics 
of the concurrence, the fohowing formulae will be very useful: 

X{A®B)^ = (detc2^) (detc2S) x^, (9a) 

(T99|(cjj^A) ® ayif) = (TA ® t(p\ay aytp) = ^ (Trc2A) x<p- (9b) 

2 Global evolution and continuous measurements 

The way to understand whether the two qubits interact or do not interact, 
directly or indirectly, is to look at the unitary dynamics of the two qubits plus 
their environment. In the Markov regime this can be done by starting from 
a quantum stochastic differential equation a la Hudson and Parthasarathy 
(HP equation) [20\ and this is also a clear way to introduce continuous 
mesurements |2H[T5j. 

As before the system space is Ji, while we take as environment space 
the symmetric Fock space % = r[L^(M;3)]; 3 is a complex Hilbert space, 
which will be only finite dimensional in the present paper. Let Ut = e"'*''^^, 
Ht = H^, denote the unitary (Hamiltonian) global evolution in % Ji. 
We suppose that the free environment evolution is @t = e~'*^°, the second 
quantization of the left shift, with its free Hamiltonian Eq. Then the global 
evolution in interaction picture with respect to is 

V{t) = e*t Ut = e'-^«* e"'*^^, t > 0, 

which, in the Markov regime, can be defined directly by a HP-equation. 



2.1 HP evolutions 



We fix a basis {Iz)}^^^ in the Hilbert space 3- Let a^(t) and ai{t) be the 
fundamental Bose field operators in r[L^(M;3)] and Az{t) = fQaz{s)ds, 
Ai{t) = jQai{s)ds, Azw{t) = al{s)auj{s)ds be the fundamental integra- 
tors of quantum stochastic calculus. 

Let us consider the HP-equation [20] for unitary operators on %f^'K 



dV{t) 



)dA,^(t)- ^ LlSz^dA^t) 



z,w£Z 



Z,W(iZ 



J2 Lz dA\{t) - \Hdt -\Y, LlLz dt 



zi^Z 



z&Z 



V{t)- (10) 
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the initial condition is V{0) = 1. By taking 



1. H, Lz, Szw G (bounded operators), Vz,it; G Z, 

2. H = H*, 

3. S £ 11(3 (8) ^) (unitary operators), where S = J2zw "S^" •S'^iui 



the solution of (|T0|1 is indeed unique and unitary. Every operator is identified 
with its natural extension to JC'^'K. 

By using the time ordered exponentials introduced by Holevo [22j, the 
solution V{t) can be represented as 



V{t) = ^ -i / Yl KzM{s)aUs) - E ^* (l 

I •'0 L ZW 



K 



S* 



aw[s) 

ZW 



ZW 



ds^, (11) 



where S = e~^^ , with a selfadjoint operator K on "K. 

Moreover, we have that Ut, defined by Ut := QtVit) for t > 0, and by 
Ut := C/J*t for t < 0, is a unitary strongly continuous group. So, we can 
interprete Ut as the evolution operator of a closed system, Qt as the free 
evolution of the fields and V{t) as the total evolution in the interaction 
picture with respect to Qf 

The interaction between J{ and % is regulated by the system operators 
H, Lz and Szw', the corresponding global Hamiltonian Ht is a very singular 
unbounded operator which could even encode the whole interaction just in 
the shape of its domain |23j . Anyway, thanks to representation (jll|) . the 
global Hamiltonian Ht has the heuristic expression 

Ht = Eo + Y^ Kzw 4(0) a^{Q) -Y,L*z [j^) ^^«-(0) 

ZW ZW 

+ E(/t) Lzal{Q)+H + Y.Ll( f 7"'!;0 (12) 
^\S-t)zw ^ V4(sin(K/2))^y.«, 



which allows to read more explicitly the interaction between the systems. 
In the special case L = we have 

HT = Eo + Y,Kzwa\z{Q)a^{Q)+H, (13) 

ZW 

while for K = i.e. S" = 1, we get 

HT = E^-\Y,L*zaz{Q)+\Y, Lz a\ (0) + H. (14) 
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As initial state let us take \e{v)){e{v)\ (8) po, where po G §{!K) is the 
initial system state and e{v) is the coherent vector in IK = r[L^(]R;3)] with 
argument v in L^(M; 3)- At the end it will be possible to take v only locally 
square integrable. 

Then, thanks to the properties of the HP-equation, the dynamics of the 
reduced system state 

Vit) := Troc {f/(t) (|e(t;))(e(t-)| O po) Uit)*} 

= Trac {Vit) (|e(t;))(e(t-)| O po) Vit)*} (15) 

is given |20^ll5j by the master equation ri{t) = C{t)[ri{t)] with Liouville 
operator 

C{t)[r] = -i[H{tlT]+Y,{L.{t)rLM* - i {L,(t)*Z,(t), r}) , (16a) 

L^{t) := L-^ +^^{Szw - 5zw)vui{t), (16b) 
w 

H{t) ■=H + ^Y1 l^iSl,, + ,5^^) + ^S,^v^{t) - h.c] . (16c) 

zw 

Of course the reduced evolution depends on the global dynamics (jlOp and 
on the environment initial state. But this correspondence is not injective 
at all, so that it is not enough to know the Liouvillian C to know the sys- 
tem/environment interaction. 



2.2 From the HP-equation to the SSE 

The fields which have already interacted with "K can be manipulated in 
various ways and then monitored continuously in time. In this way we 
avoid to further perturb the dynamics of "K, but, at the same time, as we 
indirectly acquire information on its state, the dynamics of turns out to be 
conditioned by the observed output. In the typical case of quantum optics 
the system is a photoemissive source and the output fields are mixed up by 
means of beam splitters and optical fibers and detected by photon counters 
(direct, homodyne, heterodyne detection) [T5]. In general, we identify a 
measurement in continuous time by a family of commuting selfadjoint field 
operators which can be chosen as follows. 

The manipulation of the fundamental fields is represented by a unitary, 
possibly time dependent, matrix Uiz{t), 

Uiz{t) Uiw{t) = 8zw, ^ Uiz{t)ujz{t) = Sij, 

iez z&z 
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and produces the new field operators 

Then, set dim3 = d+d' , we choose as observables the commuting selfadjoint 
operators (interaction picture) 

Bi{s) + Bj{s), Afcfe(s), i = l,...,d, k = d + l,...,d + d', s > 0. (17) 

The global evolution (|lUp . the environment initial coherent state \e{v)){e{v) 
and the observed fields (|17p together determine both the distribution of the 
output processes and the a posteriori dynamics of the system Ji, that is 
the evolution of as a function of the observed outputs; both of them de- 
pending on the system initial state pQ. As we observe a maximal family 
of compatible fields, the a posteriori evolution preserves the purity of the 
system states and thus the problem of dynamics and observation can be 
reduced to a classical linear SSE [STlfTGl: 



dip{t) = K{t)ip{t^)dt + Rj{t)ip{t^)dWj{t) 



fc=i 



(18) 



K{t) := -iHoit) -^Y.R^(trRj{t). (19) 
Ho{t) ■=H+'- {^Sl.L^ - LlS^.v,) , (20a) 



z&Z \ w&Z 



,vUt)], Jk{t) := Rd+k{t); (20b) 



the initial condition is ip{Q) = G ^) ll'^oll = 1- Equation ([18]) is a 
stochastic differential equation for a ^-valued stochastic process ip{t) in a 
filtered probability space, say (0, 3", (3"^), Q) , where Wj, are independent 
Wiener and Poisson processes, each A^^^ with rate A^. The solution ip[t) is 
taken continuous from the right and '^{t-) in the right hand side means that 
the value of the solution is taken before of the possible jump at time t. The 
solution ip{t) is a function of the initial condition -00 and of the trajectories 
of the processes Wj and A''^ up to time t. 
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Equation (fTSj) can be translated in the language of stochastic processes 
a{t) taking values among positive operators on JC. Indeed, if A{t, s) is the 
fundamental solution of Eq. (|18p . or the propagator from time s to t, taken 
po £ the stochastic process a{t) := A{t, 0)poA{t, 0)* satisfies the linear 

stochastic master equation 

d 

da{t) = C{t)[a{t^)]dt + ^{Rj{t)a{t^) + a{t^)Rj{ty)dWj{t) 

i=i 

d' 

+ E 

k=l 

where C{t) is the Liouville operator defined in Eqs. (jl6p . 

Starting from Eq. (jlSp or Eq. ()2ip one can get both the distribution of 
the outputs and the a posteriori dynamics of IK. 

Of course, the joint distribution of the compatible field observables Bj{t)+ 
Bi{t) and Afc(t) is given by the Born rule based on their joint projection val- 
ued measure and the initial system/field state. Anyway, it can be obtained 
directly from Eq. (jlSp or Eq. (|2ip. as it is the joint distribution of the pro- 
cesses Wj, Nk under the physical probability on (^7,3"^): 

PT(da;) = pTicoMdu), pt = Tr {a{t)} . (22) 

Moreover, by defining p{t) := when pt > 0, and by taking an arbi- 
trary state for p(t) when pt = 0, we obtain the so called a posteriori state, 
the conditional state to be attributed to the system, having observed the 
realization of all the processes Wj and Nk up to time t. Correspondingly, 
let us call a{t) the non normalized a posteriori state 

In particular, regarding the distribution of the outputs, by Girsanov 
theorem we can say that under the physical probability 

Wj{t) := Wj{t) - [ mj{s)ds, mj{t) := 2ReTT {Rj{t)p{t^)} , (23) 

J = 1, . . . , d, is a d-dimensional standard Wiener process, while Nk{t) is a 
counting process of stochastic intensity Pkit) = Tr { Jfc(t)* Jfc(t)/9(t_)}. 

As we observe the fields without introducing any new disturbance on "K, 
we have that its a priori state, that is the mean of its a posteriori states, 
coincides with its reduced state (fTSj) in absence of measurement: 

V{t) = Ep^ [p{t)] = Eq [a{t)] , t e [0, T] . (24) 

Thus, the continuous measurement gives an unravelling (with a physical 
interpretation) to the open dynamics (|16p . Of course, if we change the 
observed fields for a given global evolution and a given environmental initial 
state, we get a different unravelling of the same open evolution. 



Jk{t)a{t_)Jk{tT 



a{t.)) {dNk{t)-Xkdt) 



(21) 
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2.3 Interacting and non interacting subsystems 

Let us finally consider a bipartite system 'K = IKi !K2 with its environment 
% = r[L2(M;3)] and their HP-evolution We are interested in the case 
of no direct interaction between the two subsystems !Ki and 'K2, but, because 
of the common environment %, the two subsystems could have or not have 
an indirect interaction. 

If the global Hamiltonian Ht were bounded, we could say that "Ki and 
^2 do not interact directly if the global Hamiltonian is 

Ht = Hq + Hi + H2 + Hqi + Hq2 

where Hq = Hq G L('X) is the free Hamiltonian of the environment, Hi = 
HI G £(IKi) is the free Hamiltonian of "Ki, H2 = H2 G £^{'^2) is the free 
Hamiltonian of IK2, while Hqi E £(X ?fi) and H02 € £(3<; IK2) give the 
interaction, respectively, of Jii with % and of 'K2 with %. 

Analogously, dealing with HP-evolutions, we say that "Ki and 'K2 do 
not interact directly if, in the heuristic representation ()12p of the global 

Hamiltonian Ht, each one of the operators Kzw, Ylw ( 5^ ) and H + 

V / zw 

Ylzw ( -^"'^"^^ ^ ) Ly, is the sum of local operators. This property is 

V4(sin{_ft:/2)j ^ zw 

independent of the basis {\z)}z<^z chosen in 3- 

In the case L = 0, this means H = Hi + H2, with H^ = H^ € L{'Ki), 
and K = Ki + K2, with K/, = g £(3 ® J{f,). 

In the case K = 0, this means H = Hi + H2, with Hi = H^ £ LCKg), 

and each = L^P + L^, with lI^^ E L{:Kii). 

An important subcase is when the subsystems !Ki and [K2 do not have 
any kind of interaction, either direct or indirect. In other words, this means 
that each subsystem 'Ji^ has its own environment %i and that there is 
no interaction between ^1 and %i on one side and 'K2 and %2 on the 
other. Thus, we say that ^1 and 'K2 do not interact, either directly or 
indirectly, if there exists a decomposition 3 = 3i © 32i that is a decompo- 
sition % = r[L2(R;3)] = r[L2(M;3i)] r[L2(M;32)] =%i®%2 such that, 
chosen a basis {Iz)}^^^^ in 3i and a basis {\z)}z<^Z2 in 32 and consider- 
ing the heuristic representation (|12p of the global Hamiltonian in the basis 
zez-iVjZ2 in 3, each addendum is an operator on IKi ®%i or on 'K2(^%2- 
This means that K^w belongs L{'K() when both z,w G Zg, while it is null 
otherwise, that ( ) belongs to LCKf) when z G Zg, and that 

i7 + V.,„ L* ( -^-sin^ ^ J ig ^]2g g^]^ Qf local operators. This property 

V4(sin(A:/2)) 

is independent of the bases chosen in 3i and 32- 

In the case K = 0, this means H = Hi + H2, with H£ = H^ G L{'K(), 
and Lz G £(5Ci) for z € Zi, G £(^{2) for z G Z2. 
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Let us remark that the Liouvilhan ()16p is not enough to understand 
whether the subsystems !Ki and 'K2 do or do not interact. 



3 No direct or indirect interaction 

Let us start by the last case presented in the previous section, when the two 
qubits do not interact either directly or indirectly through a common bath, 
and let us study the role of a complete continuous measurement. We consider 
only the case S" = 1, so that we need to take Z = ZiD Z2, ZiCi Z2 = 0, 

L. = |^^^/ H = H,0t + t^H2. (25) 

As it will be useful in the following, from now on we give evidence to the 
tensor product structure of the various operators we need. Then, from Eqs. 
(fT6|) we get the Liouville operator C{t) = (g) 1 + 11 (g) C2{t) with 

Q{t)[T] := -imt),T] + ^ (l^tL: - \ [LlL,.r]\ , 

Hi{t) ■.= Hi + \Y, (vAt)L, - vS)K 

Recall that v is the argument in the environment initial coherent state. 
Moreover, in the case of a complete observation, we obtain the SSE (|18p 
with R^(t) and J^it) given by Eq. (Imll . K(t) = Ki(t) (g) 1 + 11 (g) K2(t), 

Let us start by considering a pure initial state = |V'o)(V'o|i so that 
a(t) = \ip{t)){ip{t)\ and pt = ||(/'(i)f , cf. Eqs. (HI]), dH]), Now the 

random a posteriori states are given by p{t) = \il}{t)){il){t)\ with il){t) = 
ip{t) / \\(p{t)\\ and the a priori states by rj{t) = Ep^, = EQ[(j(f)], see 

Sect. El 

3.1 The a posteriori concurrence 

By the definition of concurrence in the case of pure states ([5]), ([6]), we can 
introduce the random a posteriori concurrence 

Cp(t) = Ctpit) = }^'^,f^}2 (26) 



and the mean a posteriori concurrence 

Ep^ [C^(t)] = Eq |] , < t < T. (27) 
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By the definition of concurrence for mixed states ([7]) and of a priori states 
(|24p . we get that the a priori concurrence is bounded by the mean a poste- 
riori concurrence: 

By the linear SSE and Ito's formula we get the stochastic differential of 
Xip(t)j which we shah need in the fohowing, 

d d' 
dx^(t) = <t)dt + ^ £j{t)x^it) dWj{t) + [qkit)dNkit) + XkX^it) dt] , 

3=1 k=l 

(29) 

where 

d 

e{t) := TVc2 {Ki(i) +K2(t)}x^(t) + Y,{TRj{t)^{t)\ay ® ayRj{t)^{t)) , 

£j{t) := Ujz{t) (Trc2 L, + 2v,{t)^ , (30) 

zez 

1 

~k 

By writing 



Qk{t) ■■= —{TJk{t)ip{t)\ay(g)ayJk{t)ip{t)) - x.p{t)- 
At 



= 'Yf^zicri+rz, (31) 

i=l 

we get 

£j{t) = 2Y,U3z{t){rz + Vz{t)), (32) 



Trc2 {Ki(i) + K2(0} = -iTrc2 {Hi + H2} 

Let us stress that the operators Rj{t) and Jk{t) are not in general local 
operators, but sums of local operators. By this fact we cannot write in a 
more explicit form the coefficients e{t) and qk{t)- 

3.2 Only local detection operators 

As already said, in this section we are considering only local operators in 
the dynamics: every qubit has its own environment and there is no direct 
nor indirect interaction between the two qubits. Now we consider the case 
in which also the detection operators are local, that is 

Rj{t) = R^At) ®t or Rj {t) = t® R%t). (34) 
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This means that we observe separately the two environments. With this fur- 
ther assumption, the stochastic differential (j29p becomes the closed equation 



(35) 



d' 



K{t) = Trc2 {Ki(t) +K2{t)}+J2 + ^detc2R°{t), 

k=i j=i 

£j (t) = Trc2 i?° (t), dk (t) = detc2 R^^+k ) • (36) 

Equation can be explicitly solved and, by stochastic calculus, we 
get ^ 

[C^{t)] = % [C^(i)] = ^^p|~ / '^('^)'^^|' (37) 

c(t) := i^k - |4(t)|) - 2 E (I™^j(*))' - Re'^(*)- 
fc=i j=i 

The first important result is that c{t) does not depend on the initial state of 
the qubits, but only on the operators involved in the reduced dynamics and 
in the observation. This result is a slight generalization of the analogous 
one in Ref. By using ([3T]) . we get, by straightforward calculations, 



i=i zez 
^2 3 



dk{t) = ^-^^^ - E kd+k)i{tf, c{t) = J2 cj{t), (39) 
i=i jez 

3 2 

Cj{t) = 2Y,(^ehji{t)^ >0, j<d, (40) 

C3{t) = \\(^j{t)f -\dj-d{t)\+Y,\hdt)\^ j>d. (41) 

By the fact that c does not depend on the initial state of the qubits we 
can extend the result to the case of an initial mixed state po and we get 

C^(t) < Ep^ = Cp, expj- ^ c(s)ds|; (42) 

we assume always complete observation. Note that the mean a posteriori 
concurrence is non increasing. Moreover, 

r+oo 

/ c(s)ds = +00 ^ lim Ep, [Cp(t)] = 0, (43) 

JO t^+co 
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and, if c{t) = c > 0, the mean a posteriori concurrence decreases exponen- 
tially. 

For what concerns the a priori states 7]{t), when the master equation in- 
volves only local operators, one can have the phenomenon of entanglement 
sudden death (ESD) [3l[T2]. Note that no revival is possible for the concur- 
rence of rj{t) due to the bound given by the mean a posteriori concurrence 



Also the a posteriori concurrence, without the mean, can be studied. 
From the SDEs ([TH]) for ip{t) and ([55]) for Xip{t) i 'we can compute the stochas- 
tic differential of the concurrence C^(() = / IIV'C^)!!^; ™ terms of the 
new Wiener process (|23p . the final result is the closed SDE 



r 



fc=i 



14(^)1 



1 (dNkit) - iik{t) dt) - Cd+kit) dt 



(44) 



where the Cj{t) are given by Eq. ([10|) . the Cd+kif) by Eq. ([^T]) . the dkif) by 
Eq. dsn) and 



nj{t) := Re£j{t) - mj{t) = -2^ (Rehji{t)] (V'(t)|siV'(t)), 



i=l 



(45) 



j hji{t) ai^l if Rj{t) = R^{t)0l, 
\ hji (t) t(g)ai if Rj{t) = t(g) i?o (t) . 

The solution of the SDE ()44p is given by the stochastic exponential 



nj{s)dWj{s) 



CvW = C'Voexp<^ V / 

"E/ (cd+fc(s)+i4(s)i -/ifc(s))dsi J] n 

fc=l"^0 ' 0<s<tk=l 



nj(s)2 



4(s) 



ds 



(46) 



3.2.1 Diffusive case 



Here we consider the purely diffusive case {d' = 0). Now, in Eq. (|37p the 
decay intensity of the mean a posteriori concurrence is c{t) = Yl'j=i'^ji^) 
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with Cj{t) given by Eq. ()40p . while the random a posteriori concurrence 
reduces to 




with n.j(t) given by Eq. ()45p . Let us stress that neither cj nor rij depend on 
the trace of the operators R^{t). 

Note that, while the a priori states T](t) can suddenly loose any entangle- 
ment (ESD), this is a.s. impossible for the a posteriori state (with complete 
observation) . 

In the particular case of all the -Rj's selfadjoint there is decay of the 
a posteriori concurrence, but, thanks to the freedom in the choice of the 
matrix u, by a change of phase we can pass from this case to the case of all 
the -Rj's anti-selfadjoint, for which there is no decay for every initial qubit 
state {rij = cj = 0). Therefore, without changing the master equation, 
i.e. without changing the dynamical behaviour of the concurrence of the a 
priori state, one gets the complete entanglement protection by the choice of 
a phase in the detection operators. The case of all the anti-selfadjoint 
gives llv'Ci)!! = constant and the SSE describes two independent random 
unitary evolutions. 



3.2.2 Jump case 

Let us consider the purely jump case, i.e. d = and Jk{t) = J^{t) (8> 1 or 



Ji^{t) = 1 J0(t), for which we have Ep^ [C^^t)] = C^^e- /o 



. d' ..t . d' 

C^(t) = C'vo exp<^ - ^ / (cfe(s)+|4(s)|-/Xfc(s))ds ^ J| J| 

fc=l 0<s<tk=l 



dk{s) 



4(t) = detc2j{?(t) 



i=l 



d' ^ 3 

= 5^Cfc(t), Cfe(t) = - |4(i)|^ - \dk{t)\ + [^ki{t) 



k=l 



i=l 



> 0. 



One can check that Cfc(i) = if and only if Im I hkj{t) hki{t)j = 0, 



Re \^ik{t) hki{t)j = 0, i,j = 1,2,3. Again, in some cases, one can protect 
the entanglement by tuning the detection operators without changing the 
mean dynamics, for instance by changing the unitary matrix u{t). Let us 
give some examples. 
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A jump operator such as = hai + i/2 contributes [12] with 

■|2 



1 



,2\2 



note that this contribution is zero when Re hi = 0, while its maximum 

conti 

±\hi 



contribution is = |^|^ + |^|^ /4 — + |^|^ /16, reached when Re hi 



The jump operator = aa± + f3 contributes with = |a|^ /2. 
Let us consider the term 

7- (^cr_ •cr+ - ^{c^+(T_,«}^ +7+ (^cr+.cr_ - ^ {c^-(T+, .j^ (48) 

in the Liouvihe operator with 7+ > 0, 5 > 0, 7_ = + 7+. Three different 
choices of detection operators, but which give rise to the same dissipative 
term ()48p in the master equation, are: 

1. J_ = -^/t^o"- and J-f = ^/7+ 0"+, which contribute to c with 7-|_ +5/2; 

2. Ji = ^/tZJIcji, J2 = -y7+o"2, J3 = \/5cr_, which contribute to c with 
6/2; 

3. Ji = (^77:^0-+ + ^/7ro-_), J2 = (^77:^0-+ - v^c7_), which 
contribute to c with ^ (-v/T^ ~ 



Note that ^ {^/T^ — \/7+) < | ^ 7+ + f • Given the dissipative term 
in the Liouville operator, the choice (3) is the best one to slow down the 
disentanglement fn[ Eq. (19)]. 

For what concerns the random a posteriori concurrence, if C^(o) > 0, 
can vanish only if (ifc(s) = detc2 J^(s) = for some k and some s, 
as in the case of a±. In the jump case we can have ESD for some trajec- 
tories, eventually for all trajectories. The exponential decay of the mean 
concurrence is due to the randomness of the time of death. 



3.3 An example with general detection operators 

Let us now consider a concrete model of non interacting qubits plus an 
environment. We want to show, in a very simple model, how the mere 
choice of the detection operators changes the behaviour of the a posteriori 
concurrence and how much this behaviour is different from the one of the a 
priori concurrence. The Liouville operator is fixed, but different choices of 
detection operators are studied. 

As staring point ([25]) we take Z = {1, 2}, v{t) = 0, 

Li = Li®l, L2 = 1®L2, Li = -L2 = W^cJx, 7 > 0, 
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The Liouville operator turns out to be 





we can also write 



'^—[(yz®'^ + '^®cr:,,ri\--jr] 



+ ^(cJx(g)lr/cTa;(8)l + ll(8)cra;r/l(g)c7x). (49) 



The master equation with Liouville operator ()49p with loq ^ {) has a 
unique equilibrium state given by ryeq = 1/4. When wq = 0, we have more 
equilibria, the statistical operators which are diagonal in the canonical basis. 
In any case the equilibrium states are separable. 

3.3.1 Concurrence of the a priori state 

If one writes down the master equation with Liouville operator ()49p . one 
sees that it decomposes in subsystems of equations which can be solved 
analytically. However, to simplify the analysis of the dynamics and the 
computation of the concurrence, it is worthwhile to consider the subclass of 
the "X" states given in Section 11.21 By checking the master equation with 
generator (|49p in the canonical basis, one can see that the class of X states 
is preserved. 

Case Wo ^ 0. By the fact that there is a unique equilibrium state propor- 
tional to the identity, we get 



Then, if the initial X state has positive concurrence, it exists a finite time 
to > ^ for which Cp{t^-^ = and we have entanglement sudden death. 

Case (jQ = 0. In this case there is not a unique equilibrium state. As 
we shall see, the a priori concurrence is always limited by the exponential 
decay ([50|) (local detection operators, difi'usive case); one can also check 
that this limit is saturated when the initial state is a Bell state. But we 
can have also ESD; for instance, take as initial state pQ = |V'o)(V'o|) V'o = 
^ (|10) + i|01)) = i±i (|/3i) + /32)), which is again an X state. By solving 
the master equation and computing the concurrence by formulae dH) we find 
ESD at the time = — ^ In (\/2 — l); moreover, for t G [0,ii3] the a priori 

concurrence is given by C,^(j) = ^ (l + e"'''*)^ — 1. 



lim p23(i) = ^ lim pii{t) = 0, 
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3.3.2 Local detection operators 



We start by considering local detection operators. The unitary matrix u 
which fixes the observed fields in Sect. 12.21 is taken to be Ujz = ^jz^^'^\ 
4>j € [0,27r]. Then, the detection operators (|20bp . reduce to 



R 



R2 



Diffusive case. Let us start by an observation of homodyne/heterodyne 
type: d = 2, d' = 0. Then, by Eqs. (I38l)-(jl2]), ([IS]), (gT]) we get the a 
posteriori concurrence 



nj{s) ds 



and the mean a posteriori concurrence Ep^ [Cp(t)] = Cp^e where 



< c = 7 

ni{t) -- 
n2{t) -- 



[C0s4)lf + (cOS(/>2)^ 



<27, 

27 cos 01 {tiJ{t)\a^®^^lJ{t)), 
2^ cos 02 {'4j{t)\t®a^^{t)). 



The important feature of this model is that it shows the dependence on the 
measuring phases: the decay constant c can take any value in the closed 
interval [0,27]. Note that c does not depend on loq. Finally, by the bound 
for the a priori concurrence, we get 



(50) 



Jump case. Now let us consider a counting observation, with the same 
detection operators: d' = 2, d = 0, = Rk- From (jl6]) we can check 
that the a posteriori concurrence turns out to be non random and constant: 
C^(t) = Ctpo- This is due to the fact that the jump operators are proportional 
to local unitaries. Thus, any initial entanglement can be perfectly protected 
just by a proper monitoring of the environment. Let us stress that the a 
priori concurrence always vanishes for long times and sometimes even in a 
finite time. 



3.3.3 Non local detection operators 

We give now an example of detection with non local operators for the same 
non interacting qubits. Now we measure in a non local way the environments 
of the qubits, but we do not change the interaction with the environments 
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and, thus, their a priori dynamics. We consider only the diffusive case 
(d = 2, d' = 0) and we take the unitary matrix u of Section [2.21 to be 

then, we get 



R2 = (ie^'^a^ 1- ie-'-^a a^) . 

By particularizing the general formulae of Sect. 13.11 we obtain that the 
stochastic differential of Xip{t) does not contain the white noise term and we 
have 

X^{t) = -lX^{t) + 7e''^^?(t), V{t) := {T^{t)\a, ® a-Mt))- (51) 
Again by stochastic differentiation, we get 

V{t) = 7e2'^x^(t) - lV{t) - iiOo£{t), (52) 
£{t) := (T(/7(t)| {a,(^t + t^ a,) <f{t)). 

By differentiation of £ we get more complicated expressions, including terms 
with stochastic differentials. Anyway, from Eqs. ([5T]) . ([52]) we obtain 

±P(i) = e-^±*(x^, ±P(0)) Tio^o fe-^^'^'-'^£{s)ds, (53) 

Jo 

7± := 7 (l ± e''') ■ 

In this model one can have a variety of behaviours for the mean concur- 
rence, such as revivals and creation of concurrence in the long run. Let us 
see this in the simplest case. 

The case ujq = 0. In this case we have 

X^it) = ^e-^+*(xv,o +^(0)) + ^e-^-*(xv,o -^(0)). (54) 

Being non random, by Eqs. (f26|) . ([27]) . we get Ep^ [C'i/)(t)] = |X(/3(t)|; for all 
T>t. 

If e^'^ 7^ ±1, we get Re 7-1- > 0. Then, the mean a posteriori concurrence 
decays exponentially at long times, but, depending on the initial state of 
the qubits, it can have also revivals. For instance, by taking ipQ such that 
X^o ~ and 'D{0) 7^ 0, we have 
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If e = 1, we get 7+ = 27, 7_ = and 

\xM = \ |e-'^*(xVo +^(0)) + (xv-o -^(0))| • 

So, depending on the initial state of the qubits, some concurrence can survive 
(entanglement protection) or can be created in the long run (entanglement 
generation). The case e^^^ = — 1 is similar. 

4 An example with indirect interaction 

In this last section we consider the case of indirect interaction between two 
qubits and, by means of an explicit model, we show that a very extreme sce- 
nario can occur: the interaction with the environment completely destroys 
any entanglement between the qubits, if no measurement is performed, while 
the same interaction generates maximally entangled states, independently 
of the initial state of the qubits, if the environment is simply continuously 
monitored after the interaction. Indeed, while in the long run the a priori 
state of the qubits becomes maximally chaotic, and thus separable, their a 
posteriori state becomes maximally entangled for every output of the con- 
tinuous measurement. 

We consider a couple of qubits '}i = "Ki ® interacting with a sort of 
continuous flow X = r[L^(M; 3)] of quadruples of qubits 3 = 3i(8>32<8'3'i<8'32- 

Let us denote by {|i)}j=o,i the canonical basis in !Ki = = 3i = 32 = 
3'i = 32 = and then let us introduce the flip operator in 3^ ^3 

Let us choose in 3 the basis generated by the Bell bases in 3i (X" 32 and 
in 3'i <8) 32, that is {|/3^ O /3x')}x,x'=o,...,3- 

We consider the HP evolution (|lUp generated by the interaction H = 0, 

L = 0, S = FiF2=F2Fi = e-'l(^i+^2-2)^ k = ^{Fi + F2 - 2), 

where every operator is identified with its natural extension. Roughly speak- 
ing, when a quadruple of qubits 3 belonging to the continuous flow interacts 
with the couple of interest !K, the first two qubits of the quadruple 3i CS) 32 
exchange their joint state with !K, while the other two qubits 3i <^ 32 are 
simple witnesses. Then 

Sixx')iyy') = Tr3 [{\f3y (E) I3y>){/3x M ajl)^] = \f3y){l3x\ Sx'y' 

and the Hudson-Parthasaraty equation is 

= Yl {\^y)(^-\ - Sxy)v{t)dA(^,,,)^y,,){t). 

xyx' 
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Therefore, there is no direct interaction between Jfi and 'K2 as K = ■^{Fi + 
F2 — 2) with Fi involving only Jfi and F2 involving only Ji2- Let us remark 
that this is just one of those cases where the whole interaction is encoded in 
the domain of the global Hamiltonian Ht- Indeed, [23] Hj- is just an exten- 
sion of the free field Hamiltonian Eq , re-restricted to the domain of the "reg- 
ular vectors" ^ £ X "K such that axx'iO~) = J2yy' S[xx'){yy') dyy'i^^) ^ 

for all x,x' . 

For the environment we choose the initial pure coherent state \e{v)){e{v)\ 
with argument 



j— 3 

V I E 1^-) ® 1/^-) G 3 = (3i 32) ® (31 ® 32 



^(*) = a/ 7 2^1/3-)® 1/3-) g3 = (3i®32) ^ (3'i0 32), VO<t <r^ 

x=0 

where is a positive parameter and T > is our arbitrary time horizon. 
Roughly speaking, even if the qubits 3'i and 32 are not involved in the 
interaction with "Ki and 'K2, they are initially entangled with the qubits 3i 
and 32 which exchange their state with "Ki and 'K2- 

Then, if pQ is the system initial state, its reduced state at time t is 



77(t) = Tr3c \Ut\e{v)){e{v)\®poUl 
where 

Tr T7 

Lt] = v a. — ur), 

so that 

77(t)=poe-'^* + |(l-e-''*)^p fort ^00, 

and the state of !K becomes maximally chaotic and any entanglement be- 
tween Jfi and 'K2 is destroyed by the interaction with the common bath. 
The a priori concurrence goes to at least exponentially, 

C^it) < e""^ 0, for t 00, 

and, depending on the system initial state po, we can even assist to entan- 
glement sudden death. 

This can be verified by considering an X state as initial state. Indeed, 
if rjij are the matrix elements of po with respect to the computational basis 
P, we find 

C^(t) = 2max{0, Ci{t), C2{t)} , 
Ci{t) = |r/23| e-'^* - y (r/iie-* + ^(1 - e-*)) [^a^--^ + e"-*)) , 

C2{t) = |r?i4| e-'^* - J (r722e--* + ^(1 - e'-*)) (r/ase"'^* + ^(^ " e""*)) • 
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By the fact that 

hm Ci(t) = hm C2(i) = -t, 

if the initial X state po has positive concurrence, it exists a finite time to > 
for which C^{f^') = 0. The death time to can be exphcitly computed. For 
example, if po = |/3i)(/?i|, then to = 

Let us introduce now the continuous measurement. As a preliminary 
step, let us suppose we observe all the sixteen compatible processes of ob- 
servables ^(xx'){xx'){t)- Roughly speaking, we count the quadruples of kinds 
(.t.t') which have been through an interaction with the couple J{ between 
time and time t. Then the corresponding linear stochastic master equation 
for the non normalized a posteriori state a(t) is 

da{t) = £[a{t-)]dt 

+ E (^\Px'){Px\Ht-)\Px){M-^{t-)) (dNxx'it) - ^dt 



16 



in a probability space {0,,3',3^tj ^xx'{t),Q) where N^x'it), x,x' = 0, ...,3, 
are sixteen independent Poisson processes of rates A/16 under Q. 

The definitive step is to consider the measurement of the (non maximal) 
family of the four compatible processes of observables 

3 

^x'{t) = y^,^(xx')(xx'){t)- 

x-=0 

Roughly speaking, we count the quadruples of qubits, with the second couple 
of kind x', which have been through an interaction with the couple Ji be- 
tween time and time t. Then, by conditioning, we get the linear stochastic 
master equation for the non normalized a posteriori state cr(t), 

da{t) = £[a{t-)]dt 

+ E (x ('^'^(*-)) IMiPx'l - C7{t^)) (dNx'{t) - jdt 
x'=0 ^ 

in a probability space {Q,S',3't, Nx'^t),^) where Nx'{t), x' = 0, ...,3, are 
four independent Poisson processes of rates A/4 under Q. 

If N{t) = X^i^/=o ^x'{t) denotes the total counts up to time t, Tn denotes 
the arrival time of the count n and if X'^ denotes the mark of count n, the 
solution is 



a{t) 
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Then, under the physical probabihty PT(dw) = Tr {(7(r)} Q(da;), the four 
counting processes N^'it) are four independent Poisson processes of rates 
which depend on the environment initial state, and the a posteriori 
state is 

Vo, ifO<t<Ti, 



Pit) 



Roughly summarizing, a flow of quadruples of qubits 3 interacts with the 
two qubits "K. Actually only the couple 3i "S^ 32 interacts by exchanging its 
state with !K, while 3'i®32 is a simple witness which is, nevertheless, initially 
entangled with 3i®32- As a result, the couple Jf becomes entangled with the 
last couple 3'i <X) 32 with which has interacted. By counting the quadruples 
gone through an interaction with 5{ and measuring the projection valued 
measure {\/3x'){l3x'\}x'=o 3'i ® 32; we get an output with the distribution 
of a marked Poisson process and, at every count, the a posteriori state of !K 
jumps into the Bell state labelled by the corresponding mark X' . 

We can also compute the random a posteriori concurrence 



Cp(t) 



Cp„ ifo<t<ri, 
1, ift>ri, 



and we find that the a posteriori concurrence goes to 1, both almost surely 
and in the mean, 

Cp{t) ~^ Ij for t — oo, P-a.s., Vpo, 
Ep[Cp(t)] = 1 - (1 - C^(o))e-"* ^ 1, for t ^ oo, Vpo, 

while the a priori concurrence goes to 0, 

MCpit)] > C\^t) ^ 0, for t 



oo. 



Therefore, while any entanglement between the qubits is a priori de- 
stroyed by the interaction with the common bath, at the same time it is 
enough to monitor the bath in a proper way to get a maximal creation of 
the a posteriori entanglement, for any initial state of the qubits. 
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